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Abstract 

In this article, we present new results of high-order coupled cluster method (CCM) calculations, 
based on a Neel model state with spins aligned in the z-direction, for both the ground- and excited- 
state properties of the spin-half XXZ model on the linear chain, the square lattice, and the simple 
cubic lattice. In particular, the high-order CCM formalism is extended to treat the excited states 
of lattice quantum spin systems for the first time. Completely new results for the excitation 
energy gap of the spin-half XXZ model for these lattices are thus determined. These high-order 
calculations are based on a localised approximation scheme called the LSUBm scheme in which we 
retain all k-ho&y correlations defined on all possible locales of m adjacent lattice sites (k < m). 
The "raw" CCM LSUBm results are seen to provide very good results for the ground-state energy, 
sublattice magnetisation, and the value of the lowest-lying excitation energy for each of these 
systems. Bowever, in order to obtain even better results, two types of extrapolation scheme of 
the LSUBm results to the limit m — > 00 (i.e., the exact solution in the thermodynamic limit) 
are presented. The extrapolated results provide extremely accurate results for the ground- and 
excited-state properties of these systems across a wide range of values of the anisotropy parameter. 

PACS numbers: 75.10.Jm, 75.50Ee, 03.65.Ca 
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I. INTRODUCTION 



The coupled cluster method (CCM) ||, |3|, £|, [|, |6|, [7|, || has previously been applied 
@, [H| m [lj, |I§ 0, 0, [T|, 0, [21], |22j to various lattice quantum spin systems 



with much success. In particular, high-order [|TJ, [L6|, [T7|, [Lj], [2(], |21| CCM calculations using 
a localised approximation scheme called the LSUBm approximation scheme have previously 
been utilised to determine the ground-state properties of the spin-half square lattice XXZ 
model to high accuracy (see especially Refs. [|T5|, [T^). These calculations were performed 
both in a regime characterised by non-zero Neel order |], |TI], O, [TJ| in the ^-direction, and 
also in a regime characterised by non-zero Neel ordering in the xy- plane [TBI, [TP]. The results 
were found to be in excellent agreement with the known results of the XXZ model (discussed 
below). 

The spin-half XXZ model is defined by the following Hamiltonian, 

1 N 

ff =^EK + ^ + A^J , (1) 

Z i=l P 

where the index % in Eq. (HD runs over all N lattice sites, and p runs over all z nearest- 
neighbours lattice vectors with respect to i on the chain (z=2), the square lattice (z=4), 
and the simple cubic lattice (z=6). Periodic boundary conditions are also assumed. We 
note that [H, Sj] = for this model, where = YlfLi s t- The spin-half XXZ model 
contains three regimes with respect to the anisotropy parameter A. For A < — 1 there is 
a ferromagnetic regime in which the classical ground state is also the quantum-mechanical 
ground state. At A = — 1, there is a first-order phase transition to a planar regime in which 
the classical ground state is a Neel state in the x?/-plane. For the square and cubic lattices, 
the quantum-mechanical ground state in this regime is furthermore characterised by non- 
zero, in-plane, long-range ordering. For the linear chain quantum fluctuations are strong 
enough to destroy all long-range ordering and so this region in ID is sometimes referred to 
as the "critical regime." In the ID case the planar regime spans the range — 1 < A < 1, 
and at A = 1 there is another phase transition to a phase characterised by non-zero Neel 
ordering in the ^-direction for the ground state. For the linear chain, the Bethe Ansatz 



f23| , [24i |25| , |26f provides an exact solution of the quantum spin-half XXZ model, although 
for higher spatial dimensionality no exact solutions have as yet been determined. However, 
extensive approximate calculations have been carried out for this model, in particular for the 
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Heisenberg antiferromagnet (HAF) at A = 1. It is believed that the ID phase transition at 
A = 1 also occurs for higher spatial dimensionality than one. Examples of such approximate 
calculations for the spin-half XXZ and HAF models are spin- wave theory p?|, [28|] (SWT), 
exact diagonalisations of finite-sized systems |30| , exact cumulant series expansions |5T| , 
and quantum Monte Carlo |32], |33|, 133, [35[ (QMC) calculations. In particular, we note 



that such results for the spin-half isotropic HAF on the square lattice typically predict that 
approximately 61-62% of the classical Neel-ordering remains in the quantum case. 

In this article we present results of high-order coupled CCM calculations for the ground- 
and excited-state properties of the spin-half XXZ model on the linear chain, square lattice, 
and simple cubic lattice based on the systematic incorporation of many-spin correlations on 
top of a Neel model state with spins aligned in the ^-direction. In Sec. II, we describe the 
high-order CCM formalisms for both the ground and excited states. Although the high-order 
CCM formalism has been described previously [pl| , we present a very brief overview of it in 
order to give a necessary background for the new high-order formalism for the excited state. 
We also discuss in Sec. II the manner in which the CCM results for a localised approximation 
scheme called the LSUBm scheme (in which we retain all fc-body correlations defined on all 
possible locales of m adjacent lattice sites (k < m)) are extrapolated to the limit m — > oo. 
In Sec. Ill, we present the results of our high-order treatment of the XXZ model for these 
lattices, and finally our conclusions are presented in Sec. IV. 



II. THE COUPLED CLUSTER METHOD (CCM) 

In this section, we firstly describe the general ground-state CCM formalism @, |3|, |], 
H> H) Si an d then show how to apply it to the specific case of the spin-half XXZ model. 
This is then extended to deal with excited states. 



A. The Ground-State Formalism 

The exact ket and bra ground-state energy eigenvectors, |\&) and of a many-body 
system described by a Hamiltonian H, 

H\*) = E g \*)- (*\H = E g (*\, (2) 
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are parametrised within the single-reference CCM as follows: 

|M>) = e 5 |$> ; S=J2 S l°t, 
(V\ = ($\Se- s ; 5=l + X)5jCf. (3) 

The single model or reference state |$) is required to have the property of being a cyclic 
vector with respect to two well-defined Abelian subalgebras of multi-configurational creation 
operators {Cf} and their Hermitian-adjoint destruction counterparts {Cf = (Cj~)'}. Thus, 
|$) plays the role of a vacuum state with respect to a suitable set of (mutually commuting) 
many-body creation operators {Cf}, 

CT\9) = , 1^0, (4) 

with Cq = 1, the identity operator. These operators are complete in the many-body Hilbert 
(or Fock) space, 

1 = |$)($| + £C+|$)($|C7 . (5) 

The correlation operator S is decomposed entirely in terms of these creation operators {Cf}, 
which, when acting on the model state ({(7/1$)}), create multi-particle excitations on top 
of the model state. We note that although the manifest Hermiticity, (($\^ = |\l/)/(\l/|\l/)), is 
lost in these parametrisations, the intermediate normalisation condition (^|^) = ($| V I / ) = 
($|$) = 1 is explicitly imposed. The correlation coefficients {Si, Si} are regarded as being 
independent variables, even though formally we have the relation, 



(<£>\S = 





e st e 


s 




e s e s \ 





(6) 



The full set {Sj, Si} thus provides a complete description of the ground state. For instance, 
an arbitrary operator A will have a ground-state expectation value given as, 

A= = ($|5e- 5 Ae 5 |$) = A ({Si, Si}) . (7) 

We note that the exponentiated form of the ground-state CCM parametrisation of Eq. @ 
ensures the correct counting of the independent and excited correlated many-body clusters 
with respect to |<&) which are present in the exact ground state |^}. It also ensures the 
exact incorporation of the Goldstone linked-cluster theorem, which itself guarantees the 
size-extensivity of all relevant extensive physical quantities. 

4 



The determination of the correlation coefficients {<Sr,iS/} is achieved by taking appropri- 
ate projections onto the ground-state Schrodinger equations of Eq. (0). Equivalently, they 
may be determined variationally by requiring the ground-state energy expectation functional 
H({Si, Si}), defined as in Eq. (0), to be stationary with respect to variations in each of the 
(independent) variables of the full set. We thereby easily derive the following coupled set of 
equations, 

dH/dSj = ($|C7e- 5 ife 5 |$) = 0, V/ ^ ; (8) 
5H/5S! = (<5>\Se- s [H,Cf]e s \<5>) =0, V/ ^ . (9) 

Equation (|8|) also shows that the ground-state energy at the stationary point has the simple 
form 

Eg = E g ({Sj}) = ($|e- 5 #e 5 |$> . (10) 

It is important to realize that this (bi-)variational formulation does not lead to an upper 
bound for E g when the summations for S and S in Eq. ([D are truncated, due to the lack 
of exact Hermiticity when such approximations are made. However, it is clear that the 
important Hellmann-Feynman theorem is preserved in all such approximations. 

We also note that Eq. ([8]) represents a coupled set of nonlinear multinomial equations 
for the c-number correlation coefficients {Si}- The nested commutator expansion of the 
similarity-transformed Hamiltonian, 

H = e~ s He s = H+[H,S} + -[[#, S], S] + ■ ■ ■ , (11) 

together with the fact that all of the individual components of S in the sum in Eq. (|3|) 
commute with one another, imply that each element of S in Eq. (|3|) is linked directly to the 
Hamiltonian in each of the terms in Eq. ([IT]). Thus, each of the coupled equations (^) is of 
linked-cluster type. Furthermore, each of these equations is of finite length when expanded, 
since the otherwise infinite series of Eq. (|TT|) will always terminate at a finite order, provided 
(as is usually the case) only that each term in the second-quantised form of the Hamiltonian 
H contains a finite number of single-body destruction operators, defined with respect to 
the reference (vacuum) state |$). Therefore, the CCM parametrisation naturally leads to a 
workable scheme which can be efficiently implemented computationally. It is also important 
to note that at the heart of the CCM lies a similarity transformation, in contrast with the 
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unitary transformation in a standard variational formulation in which the bra state is 
simply taken as the explicit Hermitian adjoint of 

We now wish to apply the general CCM formalism outlined above to the specific case 
of the spin-half XXZ model, and we choose the Neel state, in which the spins lie along the 
z-axis, to be the model state. Furthermore, we perform a rotation of the local axes of the 
up-pointing spins by 180° about the y axis, so that spins on both sublattices may be treated 
equivalently. The (canonical) transformation is described by, 

s x -> -s x , s y -> s y , s z -> -s* . (12) 

The model state now appears mathematically to consist of purely down-pointing spins in 

4 = s t± is t 



these rotated local axes. In terms of the spin raising and lowering operators = sf ± is y k 



the Hamiltonian may be written in these local 

l N r i 

H =-7 EE 4sl P + s-s^ p + 2^sl +p . (13) 

i p L J 

In this article, we also wish to perform high-order CCM calculations for this model, and in 
order to do this we firstly define the {Cf } operators more explicitly, as Cf = sf ± sf 2 • • • sf , 
where the set-index / = {ii,t2, Each of the single-site indices is allowed to cover all 

lattice sites, although double (or greater) occupancy of any particular site in any set / is 
explicitly prohibited for the spin-half case, since (s + ) 2 = 0. We analogously define the ket- 
state correlation coefficients as Sj = <Si 1 ,* a ,...,» i for a cluster I defined above. By construction 
the coefficients 5;^...^ are completely symmetric under the interchange of any two indices. 
It is also useful to define the following two operators, 

F k = E E 1 Sk,h,—M-i <■■■ <_ x ; (i4) 

l ii— 

Gkm = E E Kl ~ 1) Sk,m,ii,~,ii-2 s t ' ' ' s t_ 2 • (15) 
I h---k-2 

The similarity transform of the Hamiltonian of Eq. (0) may now be written in terms of 
these operators. We find explicitly that H\<&) = e~ s He s \§) = {Hi + H 2 + H 3 )\$), where 



/' 



,l+p 



1 



+ 2G iji+p F i F i+p + 5*?*& P \44 + P , (16) 



H * = ^EE{f (^4 + + (G 8 , 1+p + ^^pK^p^p + ^4)} , (17) 



^3 = + G ^ P + F ^+p \ • ( 18 ) 



4^ VI 2 

i p 

We note that the similarity-transformed Hamiltonian acting on the model state |$) has now 
been written purely in terms of spin-raising operators. The problem of determining a given 
CCM ket-state equation, defined by Eq. (§) for a given cluster of index /, thus becomes an 
exercise in pattern-matching the spin-lowering operators in Cj to the terms in H contained 
in Eqs. (|I6|^T8|). This task is perfectly suited to implementation using computer algebra 
techniques, and the ensuing set of coupled, non-linear ket-state equations is then easily 
solved (for example, using the Newton- Raphson method). Once the ket-state correlation 
coefficients have been determined it is then possible to find the bra-state coefficients similarly, 
as described in Ref. [IB]]. One may finally calculate any ground-state expectation values 



that one wishes to obtain in terms of the coefficients {Si, Si}. An example of this is the 
sublattice magnetisation, given by 

2 N 

M=--<*|$>*|*> , (19) 
iV i=i 

in the rotated spin coordinates defined above. The quantity provides a measure for the 
amount of Neel ordering in the z-direction remaining in the CCM ground state, \^f), with 
respect to the perfect ordering (M = 1) of the model state. 

The CCM formalism is exact in the limit of inclusion of all possible multi-spin cluster 
correlations within S and S, although in any real application this is usually impossible to 
achieve. It is therefore necessary to utilise various approximation schemes within S and S. 
The three most commonly employed schemes have been: (1) the SUBn scheme, in which all 
correlations involving only n or fewer spins are retained, but no further restriction is made 
concerning their spatial separation on the lattice; (2) the SUBn-m sub-approximation, in 
which all SUBn correlations spanning a range of no more than m adjacent lattice sites are 
retained; and (3) the localised LSUBm scheme, in which all multi-spin correlations over 
distinct locales on the lattice defined by m or fewer contiguous sites are retained. Note that 
for this system we also make the specific and explicit restriction that the creation operators 
{C/} in S preserve the relationship that, in the original (unrotated) spin coordinates, = 
J2i s i — in order to keep the approximate CCM ground-state wave function in the correct 
(sy = 0) subspace. The number of such distinct (or fundamental) configurations for the 
ground state at a given level of approximation is labelled by Np. We denote as distinct 



configurations those which are inequivalent under the point and space group symmetries of 
the both the lattice and the Hamiltonian. 

In practice we find that the CCM equations often terminate at specific critical values 
(denoted A c ) of the anisotropy parameter which are dependent on the particular approxi- 
mation scheme chosen, such that no solution based on this model state exists for A < A c . 
At these critical points the second derivative of the ground-state energy is found to diverge, 
and the critical points have been shown to reflect the corresponding phase transition point 
in the real system. 



B. The Excited- State Formalism 



We now turn our attention to the CCM parametrisation of the excited state developed 



by Emrich []36], and as a specific example we present results for the high-order treatment 
of the excited states of the spin-half XXZ model. An excited-state wave function, |^ e )) is 
determined by linearly applying an excitation operator X e to the ket-state wave function of 
Eq. @, such that 

|* e > = X e e 5 |$) . (20) 

This equation may now be used to determine the low-lying excitation energies, such that the 
Schrodinger equation, E e \^ e ) = H\^f e ), may be combined with its ground-state counterpart 
of Eq. (0) to give the result, 

e e X e |$) = e- s [H,X e ]e s \<S>) (= R\<$>)) , (21) 

where e e = E e — E g is the excitation energy. By analogy with the ground-state formalism, 
the excited-state correlation operator is written as, 

X e = , (22) 

where the set {C/} of multi-spin creation operators may differ from those used in the 
ground-state parametrisation in Eq. (Q) if the excited state has different quantum numbers 
than the ground state. We note that Eq. fl2"2] ) implies the overlap relation ($|\l/ e ) = 0. By 
applying ($\Cj to Eq. ( |21~D we find that, 

e e Xf = mCye- s [H,X e }e s \^) ,V/^0 , (23) 



which is a generalised set of eigenvalue equations with eigenvalues e e and corresponding 
eigenvectors Xf, for each of the excited states which satisfy ($|\l/ e ) = 0. Analogously to the 
ground-state case, we define the excited-state operators 

Pk = E E l : <•-•<-, > (24) 

I ii—ii-i 

Qkm = E E *(* - 1) < " " " <_ 2 • (25) 

The state, -R|$) = e~ 5 [il, X e ]e 5 |$), may now be divided into three elements, R\<&) = 
(Ri + R 2 + i?3)|$), such that (after collecting together like terms) we find the following 
expressions, 

R-l = — E EI^C^^+P + 7)Qi,i+p) + %RiRi+pGi,i+p 
i P L ^ 

"I - Qi,i+p^i,i+p Qi,i+pRi-Ri+p RiRiRi+p\ &i &i+p > (26) 



1 

2 

i „ „ r i 



+ r,RR+p^i+p S t+p + Qi,i+pRi+p s t+p f 5 (27) 



^3 = -^EEj^M+P + ^M • (2 8 ) 



We note again that p runs over all nearest-neighbour lattice vectors and that terms which 
are equivalent under a lattice translation which connects one sublattice to the other have 



been collected together in Eqs. (p6|) - (Pq) . In order to determine explicitly the eigenvalue 
equation of Eq. ( p3|) we now pattern match the configurations in the set {Cj} to the spin- 
raising operators contained in R, analogously to the ground-state procedure. For low orders 
of approximation, this may be performed readily by hand, although for higher orders of 
approximation may one must again use computational methods. 

We note that the lowest-lying excited states for the XXZ model lie in the Sj, = +1 
and sfi = — 1 subspaces with respect to the "unrotated" ground state. We thus restrict 
the "fundamental" clusters in the set {C/} used in Eq. (^) to be those which reflect 
this property, and the number of such fundamental configurations for the excited state 
at a given level of LSUBm approximation is denoted by Np e . In order to solve fully the 
eigenvalue problem at a given value of A, we note furthermore that one must firstly fully 
determine and solve the ground ket-state equations in order to obtain numerical values for 
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the set {Si} which are then used as input to the eigenvalue problem of Eq. (|23|). The 
level of LSUBm approximation is also explicitly restricted to be the same for the ground 
and excited states, such that the calculation is kept as systematic and self-consistent as 
possible. Finally, we note that for the specific case of the spin-half XXZ model considered 
here, we only consider the lowest eigenvalue of Eq. (^T|) in order to calculate the excitation 
energy gap. We furthermore note that this eigenvalue is not specifically restricted by our 
truncation procedure to be a real number as the generalised eigenvalue problem of Eq. (|21~D 
is not constrained to be symmetric. However, in practice we find that this eigenvalue is real 
at every value of A considered, for each lattice, and at each level of LSUBm approximation. 
This provides a rather stringent check on the approximation scheme. 

C. Extrapolation of CCM Results 

We present results below of high-order CCM calculations for the spin-half XXZ model 
for the linear chain, square lattice, and simple cubic lattice. In particular, we determine the 
ground-state energy, the sublattice magnetisation, and the lowest-lying excited-state energy 
for these models using the LSUBm approximation scheme. It is clearly useful to be able to 
extrapolate the "raw" CCM results, at each value of A separately, to the limit m — > oo. 
In the absence of any rigorous scaling theory for the results of LSUBm approximations, 
two empirical approaches are outlined below and utilised. The first approach, denoted 
to as "Extrapolated(l) CCM," assumes a leading "power-law" dependence of the LSUBm 
expectation values with m. The value for the ground-state energy, sublattice magnetisation, 
or excitation energy (at a given value of A) is denoted by yi for a given value of m where 
Xi = 1/m. The index i denotes the i th data element of p such elements, although we note 
that i and m do not have to be explicitly equivalent. The leading "power-law" extrapolation 
is now described by, 

Ui = a + bx\ . 

We plot log(xj) against \og(yi — a) and the best fit of the data set, {xi, yi}, to the power-law 
dependence given above is obtained when the absolute value of the linear correlation of these 
points is maximised with respect to the variable a. This value of a is then assumed to be 
the extrapolated value for y^ in the limit m — > oo. 
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The second such extrapolation scheme of the LSUBm data, denoted the "Extrapolated(2) 
CCM" scheme, uses Pade approximants in which the data set is modeled by the ratio of two 
polynomials, given by 



Vi 



J2j=0 a j X i 



This furthermore implies that, 



a + aiXi + a 2 x ri H h a k x i 



Vi 



(hxi + b 2 x\ H V bix\)yi . 



We now wish to determine the coefficients aj and bj in order to find the polynomials in the 
equation above, and the above equation is rewritten in terms of a matrix given by, 



2 



1, Xi, X 

1, x 2 , x\ 



•, x 1 



k 



2 I \ I 

Ziyi, zfj/i, • • •, x\yi \ i 
X2V2, xlv2, ■ ■ ; x l 2 y 2 



a i 

a k 



1/2 



(29) 



Zp|/p, z;?/p, ■ • ■, x l p y p J \ ~bi J \y P J 
(Note that k + I + 1 = p.) This problem may be easily solved by inverting the matrix 
in order to determine the coefficients {a{\ and {bi} at each value of A separately. In the 
limit, m — > oo, it is seen that the extrapolated value of yi is given by ao- We furthermore 
note that the case / = corresponds also to a "least-squares" fit of the p data points to a 
(p — I) th -order polynomial. 

Previous extrapolation results [[B5| for the HAF on the square lattice fitted y^ against 



,1 



1/m 2 for the ground-state energy, E g /N, and the critical points, A c , while the sublattice 



magnetisation was fitted against X{ = 1/m. These "rules" were chosen because the CCM 
SUB2-n calculations for the square lattice HAF were found to converge as a function of n 
to their full SUB2 solution in these ways, and so analogous rules were used for the LSUBm 
data. Note that results for the HAF on the square and simple cubic lattices using these 
procedures are also quoted below. For the spin-half square lattice HAF, we find that the 
present rules are reasonable approximations for the scaling of the LSUBm results. 

Finally, we note that the LSUB2 results for and Xi generally fit in rather poorly to the 
asymptotic behaviour of yi and Xi as a function of i. As we are interested in the asymptotic 
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value of i/i, it is found that we obtain better extrapolated results by discarding the lowest 
order LSUB2 results where possible, i.e., for the linear chain and square lattice results. Hence 
for both extrapolation schemes used here, LSUBm results are used with m = {4, 6, 8, 10, 12} 
for the linear chain, with m = {4, 6, 8} for the square lattice, but with m = {2, 4, 6} for the 
simple cubic lattice. 



III. RESULTS 



The results for the ground-state energy, E g /N, of the CCM treatment of the spin-half 
linear chain XXZ model regime A > 1 are found to be in excellent agreement with exact 



results p3| , p4| , |25| , |26| over the whole of this regime, (and therefore no graphical plot of this 
data is given). Table | indicates the accuracy of the raw and extrapolated LSUBm for the 
spin-half linear chain HAF (A = 1) in comparison with exact Bethe Ansatz [[23|, |25|, |26| 
results. It is seen that all of the extrapolated results for the ground-state energy, apart 
from those of the extrapolated(2) CCM scheme with I = 3, are accurate to at least five 
decimal places. Indeed, the spin-half linear chain HAF is the "worst-case-scenario" within 
this regime for such CCM calculations. Thus, for A > 1 even better accuracy is obtained, 
and indeed the CCM results are exact at all levels of approximation in the limit A — > oo, 
in which case the model state is the exact quantum ground state. The results for the 
sublattice magnetisation, M, of this model using the CCM are shown in Fig. |l|. It is seen 
that the CCM results tend to the exact results as A increases, although they are non-zero 
and positive at A = 1. The raw CCM results may be greatly improved by extrapolation, 
and the extrapolated(l) CCM results of M — —0.0247 is in good agreement with the exact 
result that M becomes zero at A = 1. The extrapolated(2) CCM results do not perform 
quite as well as the extrapolated(l) CCM results, although they do present a very significant 
improvement on the raw LSUBm results. However, Fig. [I] shows that the results of both 
extrapolations are in good qualitative agreement with exact results, and we note that the 
spin-half HAF on the linear chain presents a particularly difficult challenge, since not only 
does all of the Neel-like long-range order inherent in the model state disappear at this point 
but also there is an infinite-order phase transition at this point (i.e., one in which the energy 
and all of its derivatives with respect to A are continuous). An indication of the efficacy 
of the CCM treatment of this model is also given for the lowest-lying excitation energy, as 
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shown in Fig. [| It is seen that the extrapolated results are in excellent agreement with 
the exact result, and are very much better than those results of linear spin- wave theory 
(LSWT) over the whole of this regime. Indeed, the manner in which the CCM results for 
the excitation energy behave as A — > 1 is in good qualitative and quantitative agreement 
with the exact results. This is in stark contrast to the results of LSWT in ID which clearly 
show completely incorrect qualitative behaviour for e e near to A = 1. Furthermore, Table 
| indicates that the extrapolated(2) CCM results for the excitation energy of the spin-half 
linear chain HAF are zero to three decimal places of accuracy. The extrapolated(l) results 
are also consistent with a gapless spectrum at A = 1. 

The various CCM results for the ground-state energy of the spin-half square lattice XXZ 
model are found to be in excellent mutual agreement over a wide range of A. This may 
be seen in Fig. |3], in which both the "raw" ground-state energy CCM LSUBm results for 
m = {2, 4, 6, 8} and extrapolated CCM results are plotted for this model. It is also seen from 
Fig. |3| that our results are in excellent agreement with those results of the QMC method 
given in Ref. [3l| for this model. Table [H] illustrates the CCM results for the special case 
of the spin-half HAF, and it is seen that CCM results are in excellent agreement (see also 
for example Refs. |K| |15|, [0|) with those results of SWT [^7|, p8"| , finite-sized calculations 
2T| , series expansions pi| , and QMC techniques [^, j3"5] . However, our results for 



the spin-half square lattice HAF perhaps indicate that the ground-state energy might be 
slightly lower than the best QMC estimate obtained to date of Sandvik J35I] which is, namely, 
E g /N = -0.669437(5). The extrapolated(l) CCM and extrapolated(2) CCM results are in 
also good agreement with the previous extrapolation [19|] of CCM LSUBm results at this 
point, which predict that E g /N = —0.66968. (We note again that the previous extrapolated 
value |TI| of Eg/N = —0.66968 was obtained by fitting to a quadratic in 1/m 2 .) 

As A approaches the HAF (A = 1), we find that the ground-state energy exponent, u, 
approaches the value v »s 2. We also note that the expansion coefficient of the I = Pade 
approximant (viz., the series in integral powers of 1/m) which becomes overwhelmingly 
dominant as A — > 1 is the coefficient which corresponds to the 1/m 2 term. Thus, the 
extrapolated(l) and extrapolated(2) results have independently reinforced what were the 
admittedly "naive" assumptions made in the "previous" extrapolation |19| in which the 
ground-state energy for the square lattice HAF is fitted to a quadratic function in 1/m 2 . 
(Note that similar results for the scaling of the ground-state energy are also observed for 
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both the linear chain and simple cubic lattice.) 

The various CCM results for the sublattice magnetisation of the spin-half square lattice 
XXZ model are similarly in very good mutual agreement over a wide range of A, as seen in 
Fig. f|. They are furthermore found to agree very well with the results of the best of other 
approximate methods for the specific case of the HAF at A = 1, as shown in Table |TJ, which 
typically yield results of M m 0.61. However, we note that the extrapolated(l) CCM result 
for the HAF of M — 0.557 given in Table || probably lies slightly too low. This is because 
the critical exponent v for the sublattice magnetisation appears to change very rapidly near 
to the HAF point for the square lattice, where, for example, we find that v — 1 at A = 1.06 
and that v = 0.490 at A = 1. It is expected however that the inclusion of the second- 
order term to the extrapolation using a leading-order power-law scaling would remedy this 
situation. We may see from Fig. [|, however, that the extrapolated(2) CCM results appear 
to work extremely well up to and including the HAF itself. (Note that qualitatively similar 
results are observed for the simple cubic lattice as for the square lattice.) 

The CCM results for the lowest-lying excitation energy of the XXZ model in the A > 1 
regime also appear to be in excellent mutual agreement, as shown in Fig. [5[ The extrap- 
olated CCM results seem to be qualitatively similar to results of LSWT, although they lie 
significantly lower and are believed to be much more accurate than the LSWT results in 
this regime. The mutual agreement of the extrapolated CCM results only ever falters very 
near to the HAF point (A = 1). Indeed, the spin-half square lattice HAF is believed to 
be gapless. The extrapolated(l) CCM results predict that the excitation energy gap disap- 
pears at A = 1.02. Furthermore, the extrapolated(2) CCM I = and I = 1 results strongly 
indicate that the spin-half square lattice isotropic HAF (A = 1) is gapless, with results of 
e e = —0.001 and e e = —0.020, respectively. Indeed, we believe that the extrapolated(2) 
CCM results provide excellent results for the excitation gap of the spin-half square lattice 
XXZ model across the whole of the regime, A > 1. 

The results for the critical points predicted by the CCM for the spin-half square lattice 
XXZ model are also presented in Table [H[ These results strongly indicate that the phase 
transition point is at (or very near to) the HAF point (A = 1), which is fully consistent 
with the behaviour that is believed to occur in this model. 

We have also determined completely new CCM results for the spin-half XXZ model on 
the simple cubic lattice. These results look qualitatively very similar to those for the square 
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lattice model and so are not plotted here. It is seen from Table |I| that the results for 
the ground- and excited-state properties of the cubic-lattice HAF are in good agreement 
with those results of LSWT. The results for ground-state energy of the spin-half cubic- 
lattice HAF are shown in Table [Til], and it is seen that the LSUB6 approximation is already 
quite near the converged limit. The extrapolated(l) and extrapolated(2) CCM results and 
the "previous extrapolated CCM" result for the ground-state energy are thus seen to be 
in mutual agreement with each other. For the sublattice magnetisation of the spin-half 
cubic-lattice HAF, the extrapolated CCM results give values of M = 0.836 ± 0.001, and 
this result is in excellent agreement with the result of LSWT of M — 0.843. By contrast, 
the extrapolated(l) CCM result appears to lie too low for the lowest-lying excitation energy 
of the simple cubic-lattice HAF. However, the extrapolated(2) CCM I = and I = 1 
results indicate that the spin-half cubic-lattice HAF is gapless with results of e e = 0.008 
and e e = —0.057, respectively, which is once again consistent with the expected behaviour 
of this model. 

The results for the critical point of the spin-half cubic-lattice XXZ model appear to 
converge to a value near to A = 1, which is again believed to be (or be near to) the correct 
result. (Note that with only two points an assumption must be made as to the scaling of 
A c with m. Hence, in analogy with SUB2-n results the LSUBm values for A c were plotted 
against 1/m 2 and linearly extrapolated, and so these results are referred to as the "previous 
extrapolated CCM" result in Table JTJ.) 



IV. CONCLUSIONS 



In this article new results have been presented for the ground-state properties of the 
spin-half XXZ model on the linear chain, the square lattice, and the simple cubic lattice. A 
localised approximation scheme was utilised within the CCM ground-state formalism, and 
this allowed high-order calculations to be very efficiently performed. The results were seen to 
be in excellent agreement with exact Bethe Ansatz calculations for the linear chain and the 
best results of other approximate methods for the square and simple cubic lattices. These 
results were extrapolated using a number of techniques which extend previous heuristic 
attempts at extrapolation. The results so obtained are of very high quality across a very 
wide range of the anisotropy parameter A, thus reinforcing our results at the specific point 
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for the Heisenberg antiferromagnet at A = 1. Indeed, it was found from these calculations 
that our previous heuristic extrapolations for the ground-state energy of the spin-half square 
lattice Heisenberg antiferromagnet did, in fact, provide very reasonable results, and that the 
naive assumptions made in this extrapolation were largely justified. Hence, for example, the 
previous "best estimate" of the ground-state energy of the spin-half square-lattice Heisenberg 
antiferromagnet of E g /N = —0.66968 was justified. Furthermore, our CCM calculations for 
spin-half Heisenberg antiferromagnet indicate that about 62% of the classical Neel ordering 
remains for the square lattice case, and that 83.6% of the classical Neel ordering remains for 
the simple cubic lattice case. 

A new formalism has also been introduced in order to treat for the first time the excited 
states of lattice quantum spin systems via the same high-order localised approximation 
scheme in a computationally efficient manner. The approximation level was kept at the same 
subapproximation level for both the ground and excited states, although the CCM ground- 
state correlations preserved the = property of the ground state and the CCM excited- 
state correlations produced a change of sf* = ±1 with respect to the CCM approximation 
to the ket ground state. Extrapolations were again attempted in the A > 1 regime and 
very good correspondence of the extrapolated CCM results for the lowest-lying excitation 
energy were observed for the linear chain in comparison with the exact Bethe Ansatz results. 
Furthermore, the results for the spin-half Heisenberg antiferromagnet, for all of the lattices 
considered here, strongly indicate that the excitation energy gaps of these models are zero. 

Future applications of the high-order formalism presented in this article will be to systems 
of complex crystallographic point- and space-group symmetries, e.g., those of the CaVO. 
Another extension will be to apply these methods to systems of higher quantum spin number, 
s. In particular, it would be of interest to apply the new high-order formalism for the excited 
state to the spin-one linear chain Heisenberg antiferromagnet which contains an excitation 
energy gap. The high-order formalism might also be applied to systems with electronic 
(rather than just spin) degrees of freedom. Finally, the application of the CCM to spin 
systems at non-zero temperature remains a future goal. 
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TABLE I: Results obtained for the spin-half HAF on the linear chain using the CCM LSUBm 
approximation with m = {2, 4, 6, 8, 10, 12}. Np denotes the number of fundamental configurations 
for the ground state, and Np e denotes the number of fundamental configurations for the excited 
state. The ground-state energy per spin, E g /N, the sublattice magnetisation, M + , and the lowest- 
lying excitation energy, e e , are shown. The LSUBm results for m = {4, 6, 8, 10, 12} are extrapolated 
using the two extrapolation schemes described in the text, and these results are compared to exact 
results from the Bethe Ansatz. 





Np 




E 9 /N 


M+ 


e e 


LSUB2 


1 


1 


-0.416667 


0.666667 


1.000000 


LSUB4 


3 


3 


-0.436270 


0.496776 


0.560399 


LSUB6 


9 


9 


-0.440024 


0.415771 


0.383459 


LSUB8 


26 


31 


-0.441366 


0.365943 


0.290251 


LSUB10 


81 


110 


-0.441995 


0.331249 


0.233098 


LSUB12 


267 


406 


-0.442340 


0.305254 


0.194577 


Extrapolated(l) CCM 






-0.443152 


-0.0247 


-0.0270 


Extrapolated(2) CCM (I = 0) 






-0.443154 


0.1351 


-0.0009 


Extrapolated(2) CCM (I = 1) 






-0.443152 


0.1043 


-0.0006 


Extrapolated(2) CCM (I = 2) 






-0.443153 


0.0938 


-0.0008 


Extrapolated(2) CCM (I = 3) 






-0.443171 


0.1128 


-0.0009 


Bethe Ansatz [23-26] 






-0.443147 


0.0 


0.0 



Figure Captions 

Figure 1: CCM results for the sublattice magnetisation of the spin-half XXZ model on 
the linear chain using the LSUBm approximation based on the z-aligned Neel model state, 
compared to exact Bethe Ansatz results. 

Figure 2: CCM results for the lowest-lying excitation energy of the spin-half XXZ model 
on the linear chain using the LSUBm approximation based on the z-aligned Neel model 
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TABLE II: Results obtained for the spin-half HAF on the square lattice using the CCM LSUBm 
approximation with m = {2, 4, 6, 8}. Np denotes the number of fundamental configurations for the 
ground state, and Np e denotes the number of fundamental configurations for the excited state. The 
ground-state energy per spin, E g /N, the sublattice magnetisation, M + , the lowest-lying excitation 
energy, e e , and the critical points, A c , are shown. The LSUBm results m = {4,6,8} are extrapo- 
lated using the two extrapolation schemes described in the text, and these results are compared to 
the best results of other approximate methods. 





N F 




E g /N 


M+ 


e e 




LSUB2 


1 


1 


-0.648331 


0.841427 


1.406674 





LSUB4 


7 


6 


-0.663664 


0.764800 


0.851867 


0.577 


LSUB6 


75 


91 


-0.667001 


0.727282 


0.609657 


0.763 


LSUB8 


1273 


2011 


-0.668174 


0.704842 


0.472748 


0.843 


Extrapolated(l) CCM 








-0.669695 


0.557 


-0.191 


1.001 


Extrapolated(2) CCM (I = 0) 






-0.669713 


0.623 


-0.001 


1.031 


Extrapolated(2) CCM (I = 1) 






-0.670619 


0.616 


-0.020 


1.044 


Previous Extrapolated CCM [19] 






-0.66968 


0.62 




0.96 


Linear SWT [27] 






-0.65795 


0.6068 


0.0 


1.0 


Second-order SWT [28] 






-0.67042 


0.6068 


0.0 


1.0 


Third-order SWT [28] 






-0.66999 


0.6138 


0.0 


1.0 


Series Expansions [31] 






-0.6693(1) 


0.614(2) 






QMC, Barnes et al. [32] 






-0.669 








QMC, Runge [34] 






-0.66934(4) 


0.615(5) 






QMC, Sandvik [35] 






-0.669437(5) 


0.6140(6) 







state, compared to results of exact Bethe Ansatz and results of linear spin-wave theory. 

Figure 3: CCM results for the ground-state energy of the spin-half XXZ model on the 
square lattice using the LSUBm approximation based on the z-aligned Neel model state, 
compared to quantum Monte Carlo (QMC) calculations [^ . 
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TABLE III: Results obtained for the spin-half HAF on the simple cubic lattice using the CCM 
LSUBm approximation with m = {2, 4, 6}. Np denotes the number of fundamental configurations 
for the ground state, and Np e denotes the number of fundamental configurations for the excited 
state. The ground-state energy per spin, E g /N, the sublattice magnetisation, M + , the lowest-lying 
excitation energy, e e , and the critical points, A c , are shown. The LSUBm results for m = {2,4,6} 
are extrapolated using the two extrapolation schemes described in the text, and these results are 
compared to linear spin- wave theory (SWT) results of Ref. [27]. The results denoted as "previous 
extrapolated CCM" (see Sec. Ill) use the 1/m 2 extrapolation "rule" for the ground-state energy 
and critical point in order to extrapolate the raw LSUBm result. However, we note that they have 
not previously been published for the spin-half cubic lattice cases. 





Np 


N Fe 


E g /N 


M+ 


e e 


A c 


LSUB2 


1 


1 


-0.890755 


0.900472 


1.873964 




LSUB4 


8 


7 


-0.900434 


0.867849 


1.117651 


0.690 


LSUB6 


181 


223 


-0.901802 


0.857192 


0.787066 


0.843 


Extrapolated(l) CCM 






-0.9026 


0.837 


-0.610 




Extrapolated(2) CCM (I = 0) 






-0.9018 


0.836 


0.008 




Extrapolated(2) CCM (I = 1) 






-0.9036 


0.836 


-0.057 




Previous Extrapolated CCM 






-0.9028 






0.965 


Linear SWT [27] 






-0.896 


0.843 


0.0 


1.0 



Figure 4: CCM results for the sublattice magnetisation of the spin-half XXZ model on the 
square lattice using the LSUBm approximation based on the z-aligned Neel model state. 

Figure 5: CCM results for the lowest-lying excitation energy of the spin-half XXZ model 
for the square lattice, using the LSUBm approximation based on the z-aligned Neel model 
state, compared to results of linear spin- wave theory. 
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